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D ' Abstract 

Several classes of gravitational backgrounds in 3 + 1 dimensions have been 
proposed as holographic duals to Lifshitz-like theories describing critical phe- 
^ \ nomena in 2 + 1 dimensions with critical exponent z > 1. We numerically ex- 

plore one such model, characterized by a temperature T and chemical potential 

T^lj- . fi, and find how to embed these solutions into AdS for a range of values of z. 

We find no phase transition going from the T <^ fj, to the T ^ /j, regimes, and 

t^ ' find that the solutions smoothly interpolate between the Lifshitz-like behav- 

ior and the relativistic AdS-like behavior. Finally, we exploit some conserved 
quantities to find a relationship between the energy density £, entropy density 
s, and number density n, if = , {Ts + fin). We show that this result is ex- 
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pected from general scaling arguments, and generalizes to i5 = -^^ {Ts + fm) 
for a theory dual to AdSrf+2 (Poincare patch) asymptotics with a local U{1) 
gauge invar iance. 



1 Introduction 

Since the Maldacena conjecture [1], AdS/CFT has become an important tool used 
to study non perturbative aspects of field theories (for a review see [2]). Tradition- 
ally, this has been used to study 3 + 1 superconformal field theories in a particle 
physics context [3], however considerably less symmetric theories can be studied us- 
ing holographic techniques [1]. It has also been fruitful to study toy models [5] that 
imitate certain features of string theory backgrounds, and seem to contain much of 
the physic^j. 

From the 3 + 1 dimensional cases, it has been learned that finite temperature 
corresponds to the presence of a horizon in the gravitational dual. Further, U{1) 
gauge symmetries in the bulk correspond to conserved number operators in the dual 
field theory. Therefore, to study a field theory at finite temperature and chemical 
potential in a holographic setup, one studies charged objects with horizons in the 
bulk: charged black holes [6]. 

Holographic techniques have recently been applied to lower dimensional non- 
relativistic systems as well [7H9] (for reviews see [lO]). In particular, much effort has 
gone into describing quantum critical behavior of these theories. Quantum critical 
systems exhibit a scaling symmetry 

t -)■ A^t, Xi -)■ Xxi (1) 

similar to the scaling invariance of pure AdS {z = 1) in the Poincare patch. From a 
holographic standpoint, this suggests the form of the spacetime metric 

ds^ = L^ (r^'de + r'^dx'dx^ij + '^\ , (2) 

where the above scaling is realized as an isometry of the metric along with r — )■ A~^r 
Q The above metric has no symmetry that mixes time and space, although there are 
also models that contain a Galilean (or other non-relativistic) symmetry [8|[T3|[T^. 
However, here we will be concerned with models that admit the metric ([2]) as a 
solution. 

Often, a good place to begin studying any system is to write down a toy or 
"phenomenological" model [71[8] to study generic properties. One may then consider 
possible embeddings into a more fundamental theory [T3 |[T5] - [T7] . such as a string 



^Of course one always prefers string theory models where one can directly state what the weak 
coupling degrees of freedom are. 

^Earlier studies of these metrics in a "brane world" scenario appear in |llj . It is also interesting 
to note that this metric, as well as several other "non-relativistic" metrics, are coset spaces [12]. 
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theory, where more information is known about the weakly coupled physics. Two 
toy models that admit the metric ([2]) have Lagrangians given by [H 



^'=16^ / ci^Xv/=^ ( i? - 2A - 1^^ - ^^2 ] ^3) 



(we call this model S' for short) and 

S = Y^ I d'x./^g {R-2A- 2(V0)2 - e'^^^JF') . (4) 

(we call this model S for short) where in either case J-" = dA is a two form field 
strength. There are also examples of actions with R^ corrections that admit solutions 
of the form (^ [2U]. 

There are several differences between the above models. First, the Lifshitz solu- 
tion to model S' is truly invariant under the Lifshitz rescaling symmetry, while the 
solution with metric (|2]) of model S has a logarithmically running dilaton, and so we 
call this a "Lifshitz-like" geometry. Note also that S has a U{1) gauge symmetry 
while S' does not, and so the field theory dual to S has a conserved particle number 
A^ with chemical potential /i. Also of note is that model S admits an exact black 
brane solution [19] that asymptotes to the metric ([2]), which can also be generalized 
to higher dimensions [21]. Finding black brane solutions for action S' has proven 
more difficult, and one often needs to resort to numeric methods [22142^ . but not al- 
ways [23] (an analogous analytic statement for an R^ extension may be found in [2B] . 
and one should also see [27] where a certain extension to this model admits an an- 
alytic black hole). For extensions and variations to the basic model S' , see [261128] , 
and for the holographically renormalized action, see [29j . 

In our current work we study the action 5*, also studied in [T9ll30p31] (for variations 
and extensions, see [^TUS^HTT] ). In the Lifshitz-like solutions, the gauge field, and 
in fact the term e^°"^T'^^—g in the Lagrangian, diverges. One immediately suspects 
that one must add a counter term on the boundary of the form * J-" A A^ and one 
would expect this to change the ensemble from a fixed f/(l) potential to a fixed f/(l) 
charge (again, fixed A^ in the boundary theory). One can also consider a fixed charge 
ensemble to simply be a limit where the chemical potential /i is large compared to 
other scales in the theory. Indeed, this is realized holographically by a geometry that 
interpolates between the Lifshitz-like solution and AdS; the large \i limit corresponds 
to a limit that leaves only the Lifshitz-like geometry. 

We now consider such a solution. We take that the full geometry is the Lifshitz- 
like geometry in the r — > limit, and is AdS for r — > oo. This type of solution 
was already numerically studied in [31], and is essentially a relativistic (actually 



conformally invariant) UV completion of the theory. In the AdS asymptotics cf) goes 
to a constant and A goes to a constant, and so the growth in both (p and A is cut 
off in this background. The size of e^'^At is essentially cut off by the scale where 
the theory becomes like AdS. Conversely, the scale in the geometry is determined 
by e^'^At, and we regard this as the chemical potential /igcom = A^-^^ (where L is the 
radius of AdS). Therefore, for such a geometry, we expect to have two regimes for 
energy scales: E <^ n and E ^ fi. For the regime E <^ fi, the quantity /i provides a 
cutoff for new physics, i.e. it is the scale at which one can excite quanta associated 
with the conserved particle number, and can be considered the scale leading to the 
RG flow of in the Lifshitz-like background. Further we can see that if we perform 
the rescaling (t, Xj,r) — )■ (A~^t, A~^Xj, A^r) for large A, we are zooming in on the 
Lifshitz-like part, and this is scaling yU — )■ oo. Combined with a time rescaling, this 
limit will eventually leave only the Lifshitz-like geometry. 

There is of course another natural energy scale that one can introduce: the tem- 
perature T. From the above, we expect there to be two regimes T -C /i and T ^ fi, 
and we would like to know how the theory interpolates between these two known 
asymptotic solutions. We were initially motivated to find a possible phase transition 
when going from one regime to the other, as fi offers a new scale. Black hole phase 
transitions are typically first order; there are two solutions for a given set of boundary 
conditions, but each solution has different energy densities. Hence, during the phase 
transition, one expects a non homogeneous phase (when considering a fixed energy 
ensemble). However, we find no such phase transition when studying the system 
numerically. A possible explanation is that we are studying the vacua of a 2 + 1 
dimensional theory. As the vacua are time independent, the space of vacua for the 
theory are defined by a 2 dimensional field theory. Such a field theory cannot have 
a phase transition that breaks the global translational and rotational symmetries 
(the inhomogeneous phase mentioned above) due to the Coleman— Mermin— Wagner 
(CMW) theorem (for a recent holographic study, see [SB])- Therefore, one might 
expect no phase transition to occur (this leaves open the possibility that for greater 
values of the spatial dimension d the theory might have different thermodynamic 
behavior). In all our results, we find a smooth, monotonic change in behavior from 
Lifshitz-like behavior to AdS-like behavior, similar to the "extremal" r/j — > case 
studied numerically in [3T| . 

We also explore the consequences of having certain conserved quantities along 
the flow, much in the same spirit as our earlier work [23l|25]. Using these conserved 
quantities, we are able to show that the energy density satisfies 

S = ^{Ts + ^n) (5) 



where s is the entropy density, and n is the number density in the dual theory. We 
argue that this relation is actually expected in AdS. Essentially, the above relation 
along with the first law of thermodynamics allows one to write a pair of first order 
linear PDEs for s{T,n) and n{T,fi). The solutions to these differential equations 
have a scaling symmetry that reproduces the expected scaling dependence of AdS. 
In fact, the scaling dependence of AdS is shown to be equivalent to relation (l5j). This 
argument generalize to 

in d spatial dimensions (of the field theory), i.e. for a theory with AdSrf+2 asymptotics 
and a local f/(l) gauge theory in the bulk. We believe these arguments to be general 
enough to apply to any theory with only two scales /x and T that has a conformally 
invariant UV fixed point. 

We present this work as follows: in section 2 we perform a reduction of the action 
S along a certain Ansatz and reproduce the equations of motion from this action. 
We find that the system reduces to 4 first order ODEs after introducing certain 
first integrals (conserved quantities), and further comment on the normalization of 
one such quantity, Q. We then perform a perturbative analysis of solutions near a 
regular horizon, around the AdS asymptotics, and around the Lifshitz-like solution, 
and display the exact AdS and Lifshitz black branes. We then discuss the setup for 
the numeric integration, and comment on how to read the chemical potential from 
the background. Finally, in section 3, we present the results of our numeric analysis, 
and the analysis that leads to equation ([6]). 



2 Analysis of the model 

2.1 Reduction 

We wish to consider dilatonic black brane solutions to the equations of motion fol- 
lowing from the action 

S = j^ [ d'x./^g {R-2A- 2{V<pf - e^"*^^) . (7) 

i07r(jr4 J 

We reduce the above action on the following Ansatz 

(h = (h(r), A = e^^''^dt (8) 



where T = dA. We will get only ordinary differential equations in what follows, and 
so we define 9 = ^. We may reduce the action to a one dimensional action |j 



S = — - / 2dtdxidx2 / drLiD (9) 

IGttGa 



with lagrangian 



It can be verified that all equations of motion associated with the action ([7]) are re- 
produced by (ITU]) as long as one uses the equation of motion for C . Here, C acts as a 
Lagrange multiplier, imposing the "zero energy" condition. Further, we note that C 
allows for generic r diffeomorphisms. We will refer to changing r coordinate as "coor- 
dinate gauge" transformations, to differentiate from the f/(l) gauge transformations 
associated with A. 

There are many conserved quantities associated with the above action. First, 
there is the conserved quantity associated with the shift symmetry [A, B, C, (p, G) — > 
{A + 26i, B — 6i,C,(j),G + 26i). This can be understood as a rescaling of the time 
coordinate, and the Xi coordinates that leaves dtdxidx2 invariant, and this then 
descends to the ID action as Noether symmetry. Further, we note that there is 
the conserved quantity associated with e*^ — )■ e*^ + const which is exactly the gauge 
symmetry associated with A. Further, we note that (A, B, G, 0, G) -t- {A, B,G,(f) + 
62, G — 062) is also a symmetry. One can view this as saying that the dilaton here 
is exactly a space dependent gauge coupling, and we may absorb a normalization of 
the gauge coupling into the definition of A. In addition, we have the Hamiltonian 
constraint as always. Now we count 4 dynamical fields, and 4 conserved quantities! 
This implies that the system is completely equivalent to a set of first order differential 
equations 

_^A+2B-C^g^y + e^+2B+^A = (11) 

^A+2B-Cq^ _ ^A+2B-Cqq _ 2g-A+2B-C+2G+2a0^^ ^ p^ ^^2) 

^a^2B-Cq^ + ae-^+2B-c+2G+2a<^^^ ^ p^ (^3) 

^-A+2B-C+G+2a4,QQ ^ Q (^4) 



^we keep track of normalization for later use 



Note, we still have not fixed a particular coordinate or f/(l) gauge. While there are 
many conserved quantities above, one can see that they do not Poisson commute, 
and so the above is not an integrable system. 

Now, let us simplify things a bit, and replace dG by Q using the last of the above 
relationships. We further take a different linear combination of the above equations, 
and find 



2e^+^^-^dAdB + e^+2i?-c^^^)2 ^ ^A-2B+c-2a^Q 
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A+2B-C ( fsi\2 



(50)2 + e^+2^+^A = (15) 

2gA+2B-c7^2a</. + a(M-95)) = Do (16) 

^a+2B-cq^ + aQe^ = Pq (17) 

^-A+2B-C+G+2c.4>QQ ^ Q (18) 

where we have redefined the integration constants AVq + 2aVo = Dq, Vq = Pq for 
convenience. We note that Pq is the only f/(l) gauge dependent quantity, i.e. it 
transforms under the map e"^ — )■ e*^ + const, while the other equations do not. In the 
above, we have used the notation Q for a non normalized charge density. 
Notice that the total contribution from the J-"^ term is 

^^2 = — ^ I 2dtdxidx2 I drd{e^)Q. (19) 

167rG4 J J 

When we compactify along an imaginary time coordinate, the integral over t simply 
becomes the inverse (unitless) temperature, and we restrict to a finite patch in xi 
and X2, so that this part of the action reads 



^^^^^j^^mi;: 



^-F^ = T7-7^2^^Q[e^]|- (20) 



when evaluated in a black hole background where r/j is the radius of the horizon. We 
note that this contribution is therefore 

^^' = 16^,^^^!—^' ]'- = le^^g^^— /X,eo., (21) 

where //geom is the potential difference between the boundary and the horizon. Be- 
cause we are using unitless t, Xi, we will have to restore units with L in this expression 
via T = y, V2 = L'^AxiAx2, to find 



L' 



1 V2 1 



This term in the action should give rise to a term 2Sjr2 = |;/igeomQ'geomV2 (see the gen- 
eral discussion in [22]) where ggeom is a properly normalized charge density. Therefore, 
we identify 

Q = g.eon.l6^G4L^ ^23) 

with ggcom properly normalized (we have chosen the charge convention that when 

/igeom is positive, so is ggeom)- 

Note that in this entire discussion, we are using e*^, e^, e^, e*" which have units 
of length, and e'^ has no units. Hence, one can see that Q must have units of length, 
and so ggeom has units of j^. This is indeed the proper normalization for a charge 
density if the potential has units of length, E ~ gV2/igeom is an energy, and so g ~ ■^. 
However, when we go to holographic variables, we will want /i to have units of energy, 
and this is given by ^^2"* (this is how we map length units in AdS to energy units 
in the field theory). The conjugate variable must also map, so we find 

qgcomL^ = n (24) 

where n is now the field theory number density. As expected this has units of L~^, 
a unitless number divided by two-volume. 

2.2 Perturbation theory at the horizon 

Near the horizon, we expect a linear zero in e^^, a simple pole in e^*", and that J-" 
and the dilaton e*^ go to constants as well. In what follows, we will fix coordinate 
transformations by changing to "entropy gauge" e"^^ = L'^r'^, such that the horizon 
radius r^ is a direct measure of the entropy density. Therefore, we expand as follows 

A = In^Lr fao(r-rft)^ + aoai(r-rfe)2 H jj (25) 

B = In(Lr) (26) 

C = ln(-fco(r-r^)"^+Ci(r-r^)5 + ...)) (27) 



r 

G = ln{ao{gh + go{r-rh) + ■■■)) (28) 

(j) = ln($^+pi(r-r^) + ■■■)• (29) 

As expected, we find a constraint on the initial conditions 

1 

Co = , ''"' (30) 



Other than this, the equations of motion simply furnish the conserved quantities 

Do = ^'^^^"° (31) 

Co 

Pq = aQaogh = (32) 

Co 

Q = *!M:. (33) 

Co 

As promised, only Pq is gauge dependent, depending explicitly on gh, the constant 
mode in e*^. 

From the expansion, we can also read the temperature 

T = ^i^. (34) 

where we have restored units with L (we will sometimes refer to a unitless temper- 
ature T = TL). Recall that in the above we have assumed that Qq is picked so that 
e^ asymptotes to Lr with coefficient 1 . Another way of packaging the same material 
is to say that oq is chosen so that e^ asymptotes to ttooLr. One would then read the 
temperature as 

T = -^. (35) 

This expression is time rescaling invariant, and so can be used in any (time rescaled) 
gauge. 

Finally, we note that the Pq = gauge has some physical interpretation at this 
point. When we consider a Euclideanized time coordinate, the choice Qh = is 
necessary to make the one form e'-'dt well defined on the "cigar" geometry. In some 
sense this makes Po = the preferred gauge. 

2.3 Perturbation theory at r = oo: AdS asymptotics 

Here we expand about r = oo, and require that the solution asymptotes to AdS. 
Our perturbative parameter will always be powers of 1/r. However, note that the 
Hamiltonian with non- vanishing charge Q plugged in will be "zero" only in the large 



r limit. Because of this, we find the set of equations 

_^A+2B-C(^g^Y + e^+2B+^A = (36) 

2gA+2B-C7 ^2d(P + a{dA - dB)) = Do (37) 

^a+2B-cq^ + aQe^ = Pq (38) 

^_A+2B-C+G+2a<^^^ = Q (39) 

easier to expand. We parameterize our expansion as 

A{r) = In(Lr) + Ai(r), B{r) = In(Lr) (40) 

C(r) = ln(L/r) + Ci(r), G'(r) = \n{g,) + Gi(r) (41) 

0(r) = ln($b) + 0i(r) (42) 

where we have taken that functions with a subscript are perturbative; they fall off 
at large r. Further, we note that in the above we have already restricted to the 
"entropy gauge" where r directly measures the horizon area. 

It is a straightforward matter to plug in the above functions and simply integrate 
the equations to find a solution. While doing so, integration constants are introduced 
that shift the boundary value of gi, and $&, and so we absorb these integration 
constants into the definition of Qb and $6. We find 

Af\ A Do + AaQgb-APo Q^ 



Q_ 



G'lW = -T^^ (45) 



Mr) = 12LV3 - 4$FZV^ ^^^) 

We note that in the above, we may remove Ab at the cost of changing the definition 
of gb by a rescaling of the time coordinate. We do so, and so we effectively set Ab = 0. 
At the horizon, this is interpreted as taking a value of oq such that the function e^ 
asymptotes to Lr with coefficient 1 (rather than any other constant). We will see 
later that this corresponds to picking an appropriate value of Dq. 

In fact one can read one interesting result directly from the above. We expect 
the "energy mode" to be determined by the -^ mode in Ai or Ci. Therefore, the 
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above gives us a very important piece of information. We first set Pq = so that the 
horizon value of e^ is (and so defines a good one-form). Comparing to the horizon 
data, we see that Dq oc Ts and Qgh oc n^, where n is the number density and /i 
is the chemical potential in the field theory. This gives us that the energy density 
E = CiTs + C2/in with Ci and C2 being constants. We will have more to say about 
this in section 3. 

There is of course the black brane solution in pure AdS space given by 



A{r) = In LrJl - (—) j , B{r) = In(Lr) (47) 

C{r) = In f ^ I G{r) = ln{g,). (48) 

<P{r) = ln(<l>fc). (49) 

This solution has Q = 0, Pq = and Dq = SL'^arf^. Note that in the euclidean 
version of the above background that e'^^''^ = is the only constant value allowed so 
that no singularity occurs at the tip of the "cigar" geometry. 

2.4 Perturbation theory about Lifshitz asymptotics 

We begin by writing down the solution for the Lifshitz type solution. The equations 
read 



A{r) = \n{aLLLr'), B{r) = In(Lr), G(r) = \n{aLgLr'^^ + 9l) 

2a0(r) = ln(r"^<l>i), C{r) = In ( — ) . (50) 



r 



Above we have used the two length scales L (appearing in B{r)) and Ll (appearing 
in C{r)). The scale L is given in terms of the cosmological constant A in the same 
way as other sections, so that the coordinate gauge choice for B{r) is the same. All 
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of the constants in the above solution are given by 

,2 {z + 2){z + l)^^ 



Ll = -^ y ^' L^ (51) 

A = -^ (52) 



2 a'^ + A 

\/z — \ a'^ 



a = I ^ z = ^— (53) 



Po = 2^t= (56) 

and in the above background, the conserved quantity Dq = 0. One may read the 
above in the same way as we read the black horizon functions: the quantity L^ is 
not an independent parameter, but depends on other constants. However, this time 
it does not depend on other "horizon" information, only on information from the 
Lagrangian a,L. Note also that we could choose a U{1) gauge where the constant 
part of At is zero, 

h = 0, (57) 

at the "horizon" and we do this henceforth, i.e. the background we are dealing with 
has (jL = and so Pq = to leading order. 

We write also the conserved quantity Q in terms of the "horizon" values of the 
fields as we did before 

Q = -^. (58) 

We now expand as before, taking 

A{r) = In(aLLir^) + Ai{r), B{r) = In(Lr), G(r) = \n{aLgLr'^^) + Gi{r) 

2a(j){r) = ln(r-^<|)i) + 2a(/.i(r), C(r) = In ( — J + Ci(r) (59) 

assuming that the perturbing functions with i subscripts are small. One may solve 
all of the equations via a single power law r" Ansatz because all terms come out 
homogeneous in powers of r. Doing so, we find solutions for the perturbed equations 
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to be 



A,{r) = Cir-i-i+^ + C^r-i-i-l - f'"^ r— ^ + A, (60) 

Ci(r) = Cir-f-i+i + C^r-f-^-i + ,,^°^^ r— ^ (61) 

2; — 1 2; — 1 



2v/^^IaL(^ + 2)L2 



r-^-2 + A 



Do(^-l) + 2Po(-^-2) ^_^_, 

1 z 1 I 7 ^ 1 



i,(r) = -Ci ^ ^ r-i-^+i - C2 , r"i'^-2 (62) 

V -2 - 1 V2 - 1 

where we have defined the useful constant 

7 = ^/{z + 2){9z + 10). (63) 



The modes above were already found in [3T] . 

A few comments are in order. First, the constant Ai appearing above may be 
removed by time rescaling, and amounts to a shift in the definition of oq . Henceforth 
we take Ai = 0. Further, the terms proportional to Do are related to the mass 
parameter given in [19lll0l|3l]. Further, the last two terms in Gi may be set to 
zero via an appropriate gauge choice, as the term r~^~^ in Gi is a pure U{1) gauge 
transformation. 

Note that if we want to add an IR (r — ;■ 0) irrelevant perturbation, the only 
mode that available is r~^~^+2 [HTj. Further, it should be noted that this term 
does not affect the values of Dq or Pq. Therefore, one can imagine fiowing from the 
Po = 0, Do = Lifshitz fixed point (with gi = 0) to the Do = 0,Po = AdS fixed 
point. 

Finally, we write down the exact black brane solution in [T9 | B0 | [3T] 




(64) 
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All constants are as before, and Vh is the location of the horizon. This background 
has 

^„.o. „„.H-^!r|3±il (65) 

V -2 - 1 
(Pq "was set to by construction, see the section about the perturbation theory 
about the horizon). This allows that all of the asymptotic functions found above 
may be removed by redefinition of Vh, or U{1) gauge, or time rescaling except for 
those multiplied by Ci and €2- 

2.5 Setup for numeric integration 

From the earlier discussions, we have found that Pq = is a convenient U{1) gauge 
choice. We will want to take a particular coordinate gauge condition as well, and 
above we have always taken 

P = In (Lr) . (66) 

This completely fixes the coordinate gauge. Further, we will find it useful to reexpress 
the differential equations in terms of "correction" functions as 

A = \n{Lr) + Ai{r) (67) 

P = In (Lr) (68) 

C = In^^") +Ci(r) (69) 

G = \n{L) + Gi{r) (70) 

0=0. (71) 

What we have essentially done here is to pull off all of the units from the functions, 
and so the resulting differential equations will be unitless if we identify the correct 
units for Pq, Dq and Q and use L to construct unitless quantities. The correct 
identification is 

(72) 
(73) 

(74) 
(75) 

where the hatted quantities are unitless. Inserting these into the above differential 
equations, we see that one has 4 first order differential equations for 3 functions. 
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Do 


= DoL 


Po-- 


= PoL' 


Q = 


--QL 



Therefore, one combination must be algebraic. We find this combination and solve 
for to find 

e'"-^ = (76) 



The remaining differential equations become 

1 e^i f 4aQe^i + Do) 

de^' - ^ ^ = (77) 

5e^^ - t:-^^ -^ -2 ^ = (78) 

2 a(e^i) r"^ 



Qar^e^^e"^^ 



(79) 



We now remark on several features of the differential equations. First off, there is 
the global symmetry already mentioned (A, 5, C, 0, G) -^ {A, B,C,(f) + §2, G — 0^2) • 
This symmetry has the effect of rescaling the charge Q in the differential equations 
by Q — ^ e^'^^Q, correctly reproducing the transformation property in equation ( 176|) . 
Hence, using this global symmetry, one may in fact set Q to be any value one wishes, 
knowing that all solutions are related to this one by the global symmetry |j. One 
can go further with this line of reasoning. Note that if we take a time rescaling, this 
has the effect of taking {A, B, C, 0, G) ^ {A + S, B, G,(I),G + S). While this is not 
a global symmetry of the reduced action, it simply rescales the action, and so maps 
solutions to solutions (this is an example of a Lie point symmetry). Further, one 
can see that only Dq and Pq rescale under this transformation. Therefore, as long as 
Dq is not zero, we may scale this to any quantity we wish. Therefore, to efficiently 



^Similar considerations were noted in [21] for extremal solutions. However, here we have 2 
quantities affected by the symmetry of AdS rescaling the coordinates (t, Xi, r) — >■ (At, Xxi, A~^r), Q 
and Dq (and Pq, but this is pure gauge). Therefore, rather than having all solutions related to one 
solution, as the authors of |31] had for the extremal solutions, we expect to have a one parameter 
family of solutions that are distinct. This is essentially getting at the fact that the thermodynamics 
should be sensitive only to ^, and this is our one parameter family. 
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parameterize the horizon data that we start with, we take the inputs Dq = Sarf^ 
(which, given Vh is always possible to set via time rescaling), and Q = 1 (which is 
always possible using the global symmetry). We note, however, that this will lead to 
fields not being normalized at infinity. Then, the rescalings used to normalize them 
are written easily in terms of the asymptotic values of the fields which is the output 
of the numeric integration. This then allows us to conclude what the "correct" value 
of Q and Dq should have been. 

Using Dq = Sar'l and Q = 1 as inputs, we numerically integrate given initial data 
Co and r^. Conveniently, Cq only has a window in which it is defined (given the black 
brane in AdS and the black brane in Lifshitz solutions) 



I < .„ < ^^^i- (80) 

Physically, when cq approaches a/^, we get the usual black brane in AdS. This should 

be a zero charge, zero chemical potential solution. As cq increases to "~~\/^5 
it approaches a Lifshitz asymptotic solution, and such a solution has a divergent 
potential and finite temperature, and so should be considered a ^ — )• oo solution. 
We always include an AdS completion, and fix the value of /i, and therefore the above 
limit should be an extremal limit. Such a solution would simply be the spacetime 
with the r — )■ limit being the Lifshitz solution, and the r — )■ oo limit being AdS 
numerically studied in [31]. These expectations will bear out in our numeric analysis. 
For completeness, we include the seed functions used to supply the initial condi- 
tions at (actually near) r = Vh, 



3 



e^^ = ao (^Vr - r/, + ai (r - r^)5 + . . . j (81) 

e^^ =co^^^ + Civ/f^^+--- (82) 
Vr -rh 

e^' = ao {go{r - r^) + ^i(r - r^f + ■ ■ ■) (83) 



^One can see this as the combined function e"'^e'^ diverging in the Lifshitz solutions. This 
combination is a global symmetry invariant and is saying that fx, as we have defined it, is becoming 
large. 
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with 

coDo 9a''cl + Qclrh{2-a^) + {a'-%)rl .^^. 

«o = T' ^1 = T^ \°^) 

_ Co (27a2c4 _ Qclr^{2 + 3«2) + rl{S + Sa^)) 
Co - Co, ci - ^3 (^85j 

(3cg - r,)r^ r3(2 + a^) - 3cgr^(2 + g^) - 9cgr,a^ + 27c^a^ ,^^, 

fi'O = ' ; 9l = ' lot) J 

CqQ ^CoThQ 

As mentioned above, after performing the numeric integration we see that e"^^ 
does not asymptote to 1. Therefore, to get a physically meaningful result, we have 
to do a time rescaling. This essentially rescales the value of Dq to be the correct 
value to give the correct asymptotic. This is the first piece of output: we get the 
physical value of Dq that should be associated with the solution, and this is used to 
determine part of the energy density. 

However, now we run into a small puzzle. The physical value of e*^^ is not 
strictly determined. We may use the global symmetry {A, B, C, 0, G) -^ {A, B, C, 0+ 
62,0 — a62) to rescale this value to any value we wish. Because this is a global 
symmetry that does not involve the metric, the stress energy tensor (and therefore 
the geometry), is not determined by this number. In fact, only global symmetry 
invariants can determine anything in the geometry (this was also noticed in [31] for 
extremal solutions). We wish to consider // as a scale in the theory, i.e. it is the scale 
at which new particles can be added/excited, and so we expect this to correspond 
to some scale in AdS, and so must be a global symmetry invariant. 

One way to fix this ambiguity is to note that the gauge coupling at infinity is the 
asymptotic value of e^""^ (which also transforms under the global symmetry). Hence, 
we would like the gauge kinetic term to go to a canonical value. This can be arranged 
by using the global symmetry to scale e'^'^'^\r=oo = 1- From the beginning we should 
have expected such a statement: only fixing both the asymptotic value of e°"^ and 
the asymptotic value of e'^^ will determine the geometry. We had traded this for e^^ 
and Q in the above discussions. Therefore, one can read the output value of /i from 
the numeric integration as 



gGlgQ(/) 



^^ = — IT 



(87) 



The factor of e^^ in the denominator accounts for the fact that e*^^ transforms under 
the time rescaling used to set the asymptotic value of e^^ to be 1, and the factor of e°"^ 
in the numerator is there because of using the global symmetry to set e"'^ = 1. Let 
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us assume that e^^ goes to 1 at infinity. Then the above simply states that what we 
have done is take some bare quantities ggcom and figcom and combined them into the 
global symmetry invariant ggcomyUgcom = (e°'^|r=ooggcom) (e""'^|r-=ooAigeom)- The second 
expression is made out of global symmetry invariants, and so can define scales in the 
geometry. These are what we use to define the physical chemical potential /igcom; and 
why we will find e°"^\r=oo = 1 a convenient representative from all solutions related 
by the global symmetry. 

As with all other formulae, the above has units restored using only L via /igcom = 
Lfi, and this defines the length scale in the geometry, and so /.i = ^ defines an 
energy scale in the field theory. One should further note that Q has changed under 
this global symmetry rescaling and so becomes an output of the numeric analysis as 

/S Qinput I 1 I 

V — e"<* lr=oo ~ e°"f'\r=oo' 

Schematically, what we have then is that (co,r/i, a; Qinput = 1, -Do,input(^h., «)) is 
input, and this gets mapped to (T, /i, Dq, Q), which are determined by the asymptotic 
values of the fields. 

We expect all relevant physical information is contained in ^, given the scaling 
invariance of AdS (see the next section). One way to scan through such values is to 
start with an r^ and scan through values of Cq in the window of allowed values to find 
a fixed value of fi. Then we may scan through r^ for different values. This process 
is fixing /i and scanning through Vh to see what the output values of T are. If T 
becomes multiple valued, then we expect a phase transition, and if it is monotonic, 
there is no phase transition. Further, we expect that as r^, <C /i we will get Lifshitz 
like behavior, and for r/^ 3> /i we will get AdS like behavior. This is essentially the 
limits T <^ yU and T ^ fj, in the field theory. 

We now turn to the results of the numeric analysis we have described above, and 
our result relating S to Ts and fin. 

3 Discussion and results 

The first point we wish to address is whether there is a phase transition. More 
concretely, we wish to determine whether there is more than one black brane solution 
given a value of fi and T. We can do this by setting (x = 2 (the value is arbitrary) 
and then seeing whether the temperature is a single or multiple valued function of r/j. 
We show log-log plotqj of the temperature T = TL as a function of r/j in figure [H We 
find that we always get monotonic behavior for a wide range of values of a. Hence 



^Maple worksheets available upon request 
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HTL) -^--^ 




12 3 

Figure 1: Above we have plotted In(TL) as a function of ln(r/i) for fixed fi = 2. The 
different graphs correspond to different values of 2;: z = 1.25 (a = 4) (red dashed 
curve), 2; = 2 (a = 2) (orange dotted curve), 2; = 5 (a = 1) (green dash-dotted 
curve), z = ^ = 8.1 (a = 0.75) (blue long-dashed curve). The solid black curve 
is a plot for the pure AdS black brane In(TL) = In(^), the asymptotic value 
of all graphs in the ln(r/j) — )■ 00 limit. The slopes of the graphs approach z as 
ln(r,j) -)■ -00. 



there is no phase transition associated with going from T <^ // to T ^ /i. Further, 
we see that we get the correct asymptotic behaviors on both ends of the graphs with 
T (X Tfi for Th ^ 2 and T oc r^ for r^ <^ 2 where z = ^4^ as before. Note that 
the scale of the changeover occurs approximately at ln(r/i) = ln(2) = ln(/i). This 
justifies the association of Lfi with a scale in AdS. 

Next, we graph the functions e^^, e'^^, e*^^, e^""^ for a = 2 after adjusting Dq and 
Q to give the correct asymptotic values of e'^S e*^^, e*^^, e^""^. We plot solutions for a 
small value of Vh and a large value of Vh in figure [2l We see no qualitative difference 
in behavior in that e"^^ , e*^^ , e*"^ , e^""^ appear to be monotonic regardless of whether 
Th <^ fi 01 Th^ fi (we omit plots for r^ ~ /i, as these look very similar too). 

Another way to analyze the solutions is to look at log-log plots of the metric 
functions and fields for a small value of rh- For this, we expect to see a section 
closely approximating the pure Lifshitz spacetime, and indeed this is what we find. 
We plot an example of this in figure [31 

Finally, we turn to the matter of the energy density. From the AdS asymptotics 
we read that the energy density is 
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Figure 2: The above plots depict the metric functions and fields as a function of 
r, the top plots for a small value of Vh = 0.4, and the bottom plots for a large 



value of r/i = 40. The plots show e^ (red dashed curve) 



dCi 



[orange dotted curve) 



gGi (^grggn dash-dotted curve) and e^°"^ (blue long-dashed curve). Plots (b) and (d) 
simply show that the asymptotics are correct with the convention that e^'^'^ is set to 
1 at the boundary. 



where we have already set Pq = 0). Therefore 

, Do -+- 4:aQgi^b 



£i 



6a 



(89) 
when 



is a unitless measure of the energy, and where we have denoted gi^i, = e ^ \r, 
1 . We plot this in figure HI 
Next, we will argue on general grounds that this expression for the energy density 



2q:0| 
c- r=oo 
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Figure 3: Above we have plotted the metric functions Ai (red dashed curve), Ci 
(orange dotted curve) and the fields Gi (green dash-dotted curve), 2a0 (blue long- 
dashed curve) as a function of ln(r) for fixed (1 = 2, and horizon radius r/j = 0.2. 
We see that all the functions have an approximately linear portion in the region 
—1 < ln(r) < 0. The slopes of these linear portions can be shown to approximately 
reproduce the Lifshitz background with z = 2,{a = 2). Further, the asymptotic 
value of Gi{r) -^ lii(2) is seen, as opposed to the asymptotic value of all other 
graphs ln(l) = 0. 
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lnf£ 




Figure 4: Above we have plotted ln(£^i) as a function of ln(r/i) for fixed (x = 2. The 
different graphs correspond to different values of 2;: 2: = 1.25 (a = 4) (red dashed 
curve), z = 2 {a = 2) (orange dotted curve), z = b {a = 1) (green dash-dotted 
curve), 2; = ^ = 8.1 (a = 0.75) (blue long-dashed curve). The solid black curve 



is a plot for the pure AdS black brane ln(£^i) = In (4-^1, the asymptotic value 

of all graphs in the ln(r/j) — )■ 00 limit. The slopes of the graphs approach as 
ln(r/i) — > —00. This is merely an indication of a finite energy held in the U{1) gauge 
field, and may be considered an extremal limit with "mass = charge" in the right 
units. One can check this by graphing ^^ and seeing that it goes to -^ (independent 



of a) for Th — > 0, however we do not display such plots. 
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should be expected. First, we recall a few facts, 



Do = L Tf^a — , Q = 

Co 4 

g, = Lf. = ^^L^ T = -^^, s = ^. (90) 

AnL Co lo7rG4 

Now, in the energy density we have a term Dq r^ r\ r^ Ts and a term Qfiffc ~ /in. 
Further, (?{, is the value of the potential at the boundary (after we fix e°"^ = 1 at 
the boundary), and g^Q is a global symmetry invariant (under shifting of G and 0). 
In deriving the above, we have assumed that e"^^ asymptotes to 1, fixing the time 
rescaling. Finally, the above observations allows us to write 

S='^{Ts + ^^n). (91) 

In fact, we can derive such a relationship using only conformal symmetry. We 
note that T and /i are energies, n and s are densities, and £ is an energy density. 
Hence, we expect that in d spatial dimensions these quantities should scale as 

(£, T, /i, s, n) -^ {X^+^S, XT, A/x, A'^s, A^n) , (92) 

and so the general functional form must be 

S = T^^'f{^), s = T'g[^), n = T'h{Pj. (93) 

Next, we assume the first law of thermodynamics 

d£ = Tds + ^idn. (94) 

We take the above differentials and equate the d^ coefficient to find 

f = 9' + j;h' (95) 



and equate the dT coefficient to find 



Using (!95i) this reads 



T\l + d)f - T'^'jrJ' = T'gd - T^T^J + fiT^-'hd - fiT^^^h'. (96) 



T^+'^Z = -^ {TT'^g + iiT'^h) (97) 
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which is identical to 

(98) 



(99) 



Further, 


the relationship f|95|) 


£= / ^ {Ts + 
d+ 1 ^ 

becomes 


fm) 








d f , 


+ j;h' + h) 


1 
= 9 




or rearranging a bit 
















rjid- 




= dT'^^^h- 
= drn. 


rjid- 


-'fih' 



(100) 

The generalization of the above argument for arbitrary number of particle species is 
straightforward, and gives 

S = -——(Ts + furii), df,^s = dTni+ [df^^Ui-d^Mj] ^. (101) 

Using this, and the Euler relation E = TS — PV + fiiNi we arrive at the familiar 
S = Pd: the stress-energy tensor is traceless. 

For the case at hand, the argument is reversible, i.e. the scaling relation can be 
obtained from a relation of the sort 

S = CiTs + C2/in. (102) 

Again, we assume the first law of thermodynamics 

d£ = Tds + fj,dn, (103) 

and write all functions in terms of the thermodynamic variables T, /i and find 

d£ = {TOts + fidTn)dT + {Tdfj,s + fid^n)dfi 

= dT{CiTs + C2i2n)dT + d^{CiTs + C2fin)diJ. (104) 

Equating differentials allows us to write two first order differential equations 

Tdrs + fxdrn = dr (CiTs + Cg^m) (105) 

Td^s + fxd^n = d^ {CiTs + Cs/in) . (106) 

Taking the integrability condition for these two equations [d^ acting on the first 
equation minus Ot acting on the second) we find that 

d^^s = drn, (107) 
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which is the same relation as the one above (llOOp . This allows us to separate the 
equations and find 

(1 - Ci)TdTS + (1 - C2)i^df,s = Cis (108) 

(1 - Ci)TdTn + (1 - C2)^Ji^^n = Can. (109) 

The above differential equations have the following solutions 

s = T^g(^^y n = f^^h^-^y (110) 

(these functions are related easily to g and h above) and the integrability condition 
O107p relates h to g via 



rpl-C2 \ 


1- 


-C\T^-ci . , 


//il-^1 


/.l-^^ J 


1- 


^ 1 y 1 


[t^-C2 



h' -T^)= ' ^ 9' ^^]- (111) 



Next, we would like to know what initial set of data specifies the thermodynamics 
in the entire (T, /i) plane. More specifically, we already know that we have a solution 
for the fi = case (this is the uncharged black brane in AdS case), and would like 
to know if this, or any additional information, can specify the thermodynamics off 
of this curve. To do so, we would like to calculate the characteristic equation for 
paths in T,fi space, parametrically given by T{a), ^{a). The characteristic curves 
are those that satisfy 

(1 - COT^ - (1 - C2)/i^ = 0. (112) 

(we note that the characteristic equations are the same for both of our linear PDEs). 
The most general solution to these equations, up to reparameterization of a, are 

These curves are given by T'^^~'^2^yU~'^^~*^^) = constant. These are either hyperbola- 
like curves that do not intersect T = or // = (the degenerate case being the 
"hyperbola-like" curve T'^^~*^2)^-(i-Ci) _ g^^ ^j. i^j^ey are power law curves that go 
through the origin at T = 0,yU = (degenerate cases given by T = or yU = 0), 
depending on the particular values of Ci and C2. For us we have Ci = C2 so that 
lines of constant — are related, which matches the scaling argument above. 

Note that the curve /i = is a characteristic curve, simply setting /xq = 0. This 
means that the initial conditions provided by the curve /i = are not sufficient 
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to determine the functions s and n off of the curve // = (see for example [41]). 
Further, no amount of additional perturbative information will yield a complete 
amount of information off a characteristic curve. This is realized by the presence 
of the general function g in the general solutions flllOp . This basically comes down 
to saying something that we already knew: all relevant information differentiating 
models depends on ^ (again, # = 1 for AdS). This is reflected in the geometry by 
saying that fi defines a scale, and T defines a scale and all solutions related by rescal- 
ing (/i, T) — !> (A/i, AT) are related by the diffeomorphism (r, t, Xi) — )■ (Ar, \^H, A~^Xj) 
that leave the AdS asymptotics alone, but shift the relevant scales in the bulk. What 
we have learned is that no amount of perturbative information around a solution can 
fix the thermodynamics. This is good because there are many different holographic 
models with a U{1) gauge invariance, all of which presumably have different thermal 
behavior. 

Finally note that along a characteristic curve, the functions g and h are con- 
stant. This allows us to quickly deduce that along a characteristic curve, a — )■ aX, 
(To, /io, So, no) ^ (Ai-^iTo,Ai-^^/xo,A^iso,A^^no). For us Ci = C2 = ^. Calling 

A = A'^'^^ we find that the scaling becomes (To,/io, So,no) — ;■ (ATq, A/io, A'^So, A'^rio)- 
This is just the scaling in AdSrf+2 noted above (the scaling of S follows from (155]) ). 

To summarize, we believe that the above arguments are general enough to apply 
to a theory where the only two scales are the chemical potential /i and the temper- 
ature T, and where the theory has a conformally invariant UV fixed point. Under 
these conditions, we expect fl98|) to hold, and that the different models are param- 
eterized by the functions (say) / (|;). This information is beyond relation fl^ . i.e. 
it is information off of the characteristic curves that it defines. 

We end by mentioning some possible future directions: 

1. It would be interesting to explore the space of models with different functions 

/(#)■ 

2. One could study field theories on a. d dimensional sphere, i.e. the black hole 
geometries with spherical horizon. In these cases we expect phase transitions. 

3. Other dimensions: one could expect there to be different thermal behavior, 
possibly phase transitions, for general d, as discussed at the end of section 1. 

4. One could calculate certain correlation functions in the dual theory, and see 
how they interpolate between non-relativistic and relativistic behavior. 

5. One could also study the dyonically charged versions of these black holes. There 
should be constant dilaton solutions where the magnetic and electric charge are 
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equal (see [42] for a recent study). 

6. One could also study theories that asymptote to the Lifshitz-like background 
and are charged under a different U{1) field (essentially adding a scale in the 
Lifshitz background). However, note the curious result in [13] where an analytic 
charged solution is found, but the potential diverges. It would be interesting to 
see if there are solutions with a finite (chemical) potential. Using similar scaling 
arguments, we expect that a relationship of the type (I102p can be written in 
such a case with Ci = C2 = -^- However, we would like to back this up with 
an explicit calculation. 
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